One-loop corrections to the curvature perturbation 
from inflation 



David Seery 

Astronomy Unit, School of Mathematical Sciences 
Queen Mary, University of London 
Mile End Road, London El 4NS 
United Kingdom 

E-mail: D . SeerySqmul .ac.uk and dj s61ScEmi .ac.uk 



Submitted to: JCAP 

PACS numbers: 98.80.-k, 98.80.Cq, ll.lO.Hi 

Abstract. An estimate of the one-loop correction to the power spectrum of the 
primordial curvature perturbation is given, assuming it is generated during a phase 
of single-field, slow-roll inflation. The loop correction splits into two parts, which 
can be calculated separately: a purely quantum-mechanical contribution which is 
generated from the interference among quantized field modes around the time when 
they cross the horizon, and a classical contribution which comes from integrating the 
effect of field modes which have already passed far beyond the horizon. The loop 
correction contains logarithms which may invalidate the use of naive perturbation 
theory for cosmic microwave background (CMB) predictions when the scale associated 
with the CMB is exponentially different from the scale at which the fundamental theory 
governing inflation is formulated. 
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1. Introduction 

For some years, the theory of inflation [H El |3l HJ El El [7] has represented our most 
successful approach to the very early universe. Although we do not know with certainty 
what physical conditions applied in the very distant past, it is believed that a number 
of fundamental properties of our universe were determined at that time. Among the 
most important of these is the ensemble of small density perturbations which later 
condensed into galaxies and gave rise to the network of cosmic structure we see today. 
Therefore, one question which must be answered by any theory which purports to 
describe the very early universe concerns the origin and nature of these perturbations. 
In inflationary models, the answer is provided by use of an era of quasi-exponential 
expansion to generate small fluctuations from the vacuum. Remarkably, the properties 
of these perturbations can be calculated by applying a minimal extension of quantum 
field theory to the expanding universe [8l[9l[T0lllIl[l2l[T3l|ll|15]. 

The inflationary perturbations are generated by quantum-mechanical effects among 
field modes of wavelength k around the time when such modes are of order the Hubble 
scale, that is, k ~ aH. If the inflationary expansion is still ongoing when these 
perturbations are generated then the local Hubble scale continues to contract, so that 
eventually they are far outside the horizon, giving k -C aH. In this regime the field 
behaves in an approximately classical fashion [16]. Therefore the evolution of the 
curvature perturbation can be calculated reliably until the end of inflation, when it is 
generally assumed to settle down to a time-independent final value, unless isocurvature 
modes exist which drive its evolution on large scales [T71 [181 HH]- In this simple picture 
we have a robust prediction of the properties of the perturbations regardless of the 
details of the inflationary era. 

The above account neglects at least one possible source of trouble, namely the 
question of quantum mechanical corrections. In any quantum theory, including inflation, 
we are obliged to study the role of so-called loop corrections, which may modify or 
destroy the leading semiclassical prediction. Loop corrections are systematic corrections 
in powers of a coupling constant, which is small if perturbation theory is to be a 
good approximation. However, perturbation series are generally asymptotic rather than 
convergent, and if the coefficients of higher-order terms are large then such corrections 
need not negligible even if the coupling constant is small. In such circumstances, these 
corrections must be taken into account when performing accurate comparisons of theory 
with observational data; for example, their use is routine when interpreting the outcome 
of particle physics experiments. Modern cosmological datasets have become remarkably 
precise [20l EU [22], and observational improvements imply that future datasets will 
possess even greater constraining power. It is important that our theoretical predictions 
keep pace with observational developments. For this reason, as in particle physics, we 
should ask whether loop corrections can modify the predictions we use to compare with 
experiment. 

If such corrections are to play a role, what form should we expect them to 
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take? Models of inflation make predictions concerning the properties of the primordial 
curvature perturbation, (. This perturbation is communicated to the cosmic microwave 
background (CMB) at decoupling. The statistical properties of ( are predicted by 
supposing that the universe is geometrically close to a patch of the classical, eternal de 
Sitter universe, in which the expansion is exactly exponential. The de Sitter universe 
has a special simplicity as a solution of the Einstein equations since it is maximally 
symmetric, implying that its curvature is constant. However, as an arena for physics 
the de Sitter universe is considerably less simple [23]. A horizon surrounds any freely 
falling observer and screens a portion of the spacetime from view. This horizon is 
associated with a finite entropy, Sds, which can be thought of as representing a finite 
state space of size ~ exp S^s available to each inertial observer. One therefore expects 
that perturbative predictions made over a volume V > exp S^s (in Planck units) are 
untrustworthy, since the volume V carries more perturbative degrees of freedom than 
can actually be available in the full quantum gravity [2U [25]. An inflationary epoch 
lasting for N e- folds causes the universe to expand by a factor e^, so we expect that 
the inflationary era cannot be described in perturbation theory for an arbitrarily long 
time; in order not to violate the volume bound, it is necessary that < ^ds for reliable 
calculations. Related limits on A^ of various sorts have been reported in the literature 
[261 [271 [28]. 

The above intuitive argument implies that we should expect the coefficients in the 
loop expansion to become large when N ^ 1. For sufficiently large A^, we will have a 
breakdown of perturbation theory and thus an apparent loss of predictivity from the 
early universe. The purpose of the present paper is to explore the details of this loss of 
predictivity, quantified as a loop correction to the observable perturbation (. Although a 
proper understanding of the breakdown of perturbation theory is important as a matter 
of principle, it is also important in practice. The apparent success of the inflationary 
prediction means that it has become routine for experiments which measure temperature 
anisotropics in the cosmic microwave background to report their flndings in terms of 
constraints on inflationary models. These constraints are usually derived using tree-level 
predictions for ( which ignore the possibility of a breakdown in perturbation theory. 

The important ramiflcations which attach to a failure of the perturbative description 
mean that loop corrections have received attention from many authors, beginning with 
't Hooft & Veltman [29]. The possibility of an infrared perturbative breakdown was 
explored by Sasaki, Suzuki, Yamamoto & Yokoyama [301, [3l], and later authors computed 
loop corrections in a variety of different settings; for example, Refs. [321 |33l (Ml [351 (Ml 
[aZllMllanillOllillliailHlllillSlllelll^ We win derive an estimate 

for the loop correction for an arbitrary potential, using the slow-roll approximation to 
control the calculation. A central feature of this analysis is the use of the non-linear 6N 
expansion [511 ESI ESI EZ] ■ This expansion provides a systematic method for computing 
the evolution of ( (and its correlation functions) in terms of the fleld perturbations at 
earlier times (and their correlation functions). Moreover, as will be described below, 
the 6N formulation has many useful properties for the purposes of loop calculations. In 
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order to apply this method consistently, it is necessary to account for the presence of 
loop corrections in the field perturbations at early times. The requisite loop correction 
has recently been calculated [58], and in this paper all these corrections are assembled 
into a loop correction for (, which allows an estimate of the onset of any breakdown in 
perturbation theory in a number of popular models of inflation. 

This paper is organized as follows. In ^ the non-linear 6N formula is briefly 
reviewed, and in §3] it is used to derive the complete one-loop correction to the power 
spectrum of (. Our result agrees with that of Byrnes et ai, who recently gave an 
expression for the loop correction valid to two loops. A feature of the SN expansion 
is that one is free to construct the initial hypersurface at any time after the relevant 
k-modes have left the horizon. This freedom is discussed in §3.31 In use the formal 
5N expression to compute an estimate for the loop correction. The loop correction is 
shown to be divergent on large scales in §4.11 and must be regularized by restricting 
the computation to some finite box of size £. In § §4.2fH73] estimates for the regularized 
loop correction are given in two approximations. In §4.21 the spectrum is assumed to 
be defined by a constant tilt, whereas in §4.31 we specialise to the case of a monomial 
potential and compute an explicit, but mo del- dependent, estimate. Finally, we conclude 
with a discussion of these results in ^ 

We use natural units in which the speed of light c, Planck's constant h and the 
Planck mass Mp = (SttG)"^/^ are set equal to unity. The metric signature is (— , +, +, +) 
and the unperturbed background metric is taken to be of Friedmann-Robertson- Walker 
form 

ds^ = -dt^ + a(t)2dx ■ dx, (1) 

where t is cosmic time and a is the scale factor of the universe. It is often more convenient 
to use a conformal time variable, t], defined by r/ = J^dt/a{t). The Hubble parameter 
satisfies H = a/a, where an overdot denotes a derivative with respect to t, and inflation 
occurs whenever e = —H/H"^ < 1. 

Our theory will consist of Einstein gravity coupled to a scalar field with potential 
V{(f)), which can be arbitrary except that it must support an epoch of inflation for some 
values of the field. The background is taken to be homogeneous with spatially-dependent 
perturbations 6(j) which satisfy the smallness condition \6(j)\ <^ \(f)\. When e <^ 1, the 
field rolls only a short distance in a Hubble time, and therefore this is known as 
the slow-roll regime. When slow-roll applies, it is often a good formal approximation 
to compute in powers of e or related small quantities obtained by taking dimensionless 
time derivatives of e, although a key theme of this paper will be that such expansions are 
usually asymptotic rather than convergent. The most important related parameter is rj, 
not to be confused with the conformal time, and defined hj t] = H/H<p. This measures 
the local curvature in the potential V, and in typical models it is small whenever e is, 
although this need not be true in general. 
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2. The 6N formula 

Our theories which govern cosmological evolution at high energies are typically written 
in terms of a number of microscopic fields whose quanta are the particle species which 
populate the universe. The interesting cosmological observables are expectation values 
of products of these fields, often in the special case where the fields which participate 
in the expectation value are evaluated at equal times but at distinct spatial positions 
{xi,...,x„}. It is usually more convenient to translate to Fourier space, and work 
instead with expectation values of fields evaluated at distinct wavenumbers {ki, . . . , k„}. 

The power spectrum and bispectrum. The lowest-order expectation value is the one- 
point function, 

(50(k)) = (27r)3(5(k)0, (2) 

where O is a dimensionless quantity which depends on the time of evaluation. At tree- 
level the one-point function is always zero if we have chosen to work in a perturbatively 
stable vacuum, but it is possible that a non-zero O may be generated radiatively. It 
was shown in Ref. [58] that to one-loop order O is given by a renormalization-scheme 
dependent number, which can be absorbed into a redefinition of the background field. 
The next-order expectation value is the two-point function, 

(50(kO50(k2)) = (27r)35(ki + k2)Pih), (3) 

where P{k) again depends on the time of evaluation and is known as the power spectrum. 
It is sometimes more convenient to work in terms of the "dimensionless" power spectrum 
V, which is related to P by the rule V = k^P/2Tr'^. 
The three-point expectation value can be written 

(50(ki)50(k2)50(k3)) = (27r)35(ki + k2 + k,)B{h, k2, fcs), (4) 

where B is referred to as the bispectrum, since in typical cases it is proportional to a sum 
of terms quadratic in P [59| 160] . To express B in terms of a dimensionless quantity, it is 
conventional to introduce the momentum-dependent parameter /nl [60], analogous to 
V and occasionally called the "reduced bispectrum," which is defined to satisfy [611 [STjfl 

B{k,MM) = ~hi.{p{k^)P{k2) + P{k2)P{k^) + P{k^)P{k^)]. (5) 

Note that B (and /nl) can always be written as a function of the three scalars 
{fci,fc2,fc3}, in terms of which it does not depend directly on the relative orientation 
of the {kj}; alternatively, one can use the momentum conservation condition kj = 
to rewrite 5 as a function of only two of the momenta, but in this case it will explicitly 
depend on their relative orientation, making B (and /nl) always a function of three 
independent degrees of freedom. 

t Sometimes /nl is defined by writing C, = C,g — |/nl * Cg (in coordinate space), where C is the full 
curvature perturbation, C,g is a gaussian random field, denotes a convolution and /nl is to be thought 
of as parametrizing a possible x^-typc contribution to C,. To linear order in /nl, this expansion for C, 
in terms of C,g produces the formula ([5]). However, since there is no motivation for such a specific type 
of non-linearity in C, we follow Lyth & Rodriguez j57] in adopting Eq. ([5]) for /nl, whatever its source. 
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Tilt. The dimensionless quantities V and /nl are generally scale-dependent. Their 
scale dependence can locally be characterized by two tilt parameters n and tt-nl |62j . 
defined by 

1 dlnP dln/NL 
n-l= and tinl = , , , (6) 

d m A; dink 

where in the definition of tt-nl we take /nl to be evaluated at equilateral momenta 
ki = k, and by convention n — 1 is chosen as the tilt of V rather than These 
numbers are typically small [63]. For a limited range of wavenumbers it may be an 
acceptable approximation to take V and /nl as almost flat over the range, whereas if it 
is necessary to characterize V or /nl over a larger range then one can write 

/ k \ (""^^I'^o / k \ "nlI*o 

P«PU,y and /NC«/.cL,.y (7) 

where k^ is some pivot wavenumber which is characteristic of the range of wavenumbers 
in question, and '1^^,' denotes evaluation at this point. However, both these 
approximations break down over an exponentially large range of wavenumbers since 
in general none of V, /nl, n or tt-nl are constant. Over such an exponentially large 
range Eq. substantially overpredicts the spectrum, which typically grows only like a 
power of In k rather than a power of k [HI [M] . Where possible we will use Eqs. (l6])-([7]) 
to give mo del- independent estimates for tilted spectra, but it is necessary to make more 
detailed model-dependent estimates where an exponentially large range of scales applies 
(see 



The 6N formula. The microscopic flelds are not observable by themselves; only a 
speciflc combination of them can be observed as the gauge-invariant perturbation ( 
[65] which is communicated to the CMB. An extremely powerful and computationally 
convenient way to obtain this combination is to employ the so-called 6N formula 
[5H 155] [66| [56] |57] . Consider the Arnowitt-Deser-Misner [67] line element 

ds^ = -N dt"^ + hij (dx* + A^^ dt) (dx^' + A^^' dt) , (8) 

in which the perturbations are parametrized by a lapse function A^, the shift vector A^* 
and the spatial metric hij. For any choice of t, we can write the spatial metric on the 
three-dimensional slices of constant time in the form hij = a{ty exp{2i/j{t,x.)}6ij^ The 
quantity ip is known as the curvature perturbation associated with this slicing. 

I Note that rtNL clearly does not capture the most general variation in /nl under a change of the 
momenta, since /nl is a function of three degrees of freedom. More generally, one could define the tilt 
of /nl as a vector quantity obtained by taking the gradient with respect to the ki. 
§ This is not strictly true, because this is not the most general perturbation of a spatial slice. The 
curvature perturbation -0 amounts to a fluctuation in the trace of hij, which describes the volume 
expansion in a local region of the slice. In principle there are further vector and tensor modes 
which describe vorticity and gravitational waves. These modes are set to zero. In general this is 
not a consistent procedure, because scalar, vector and tensor perturbations mix beyond linear order. 
Therefore the presence of scalar perturbations will induce perturbations in the vector and tensor modes. 
However, we expect that these effects will not destroy the leading order prediction for C, which we can 
obtain by focusing simply on -0. 
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Consider any two hypersurfaces at times (the initial hypersurface) and t^ (the 
final hypersurface), where tc > t*. In the perturbed universe the number of e- foldings 
of inflation between these slices will vary from place to place, and satisfies 

N{t, t„x) = No{t, -> tc) + ^/'c(x) - M^), (9) 

where A'^o represents the number of e-folds which would have taken place in the 
unperturbed background. By choosing the initial slice to be fiat, so that ip^^ = 0, and 
the final slice to be uniform density, so that ipc = Cc one obtains the simple expression 
( = SN, where 6N = N{x.) — Nq is the shift in e-folds between t* and tc generated 
by the perturbations. We expect widely separated Hubble volumes to evolve locally 
like independent universes, so on large scales this shift can be computed merely by 
accounting for the variation in initial conditions from place to place at time t^, 

c(x) = j2{5M^)r (^)"iv(t„t,), (10) 

n=l 

where the 50* are the field perturbations at t*, evaluated on the flat slicing. (Note that 
(c defined by this expression has no dependence on t*.) In principle N depends on all 
the initial conditions on the initial slice, and because (p is governed by a second-order 
differential equation it will generally be necessary to specify the pair {0*, 0*} at t*. We 
will make the simplifying assumption that the slow-roll approximation applies at the 
initial time. This relates 0* to 0*, and hence it is unnecessary to include 0* separately. 

Expectation values. Eq. (fTOj) can be used to compute n-point correlation functions of 
the curvature perturbation (, by first constructing a product of n copies of (JTOl) and 
taking expectation values (■ ■ ■) in the appropriate vacuum state. After application of 
this process, and working to lowest non-trivial order in 50, it follows that 

(C(X0C(X2))C = {^^^^^} V0(X1)50(X2))*, (11) 

where the subscript c denotes evaluation on the final slice at time tc- Eq. ( iTTl) shows 
that we can compute the power spectrum at any late time t^ using our knowledge of 
correlators only at some fixed early time t*. By repeating the above process for higher 
n-point expectation values and to higher orders in 50, it is obvious that we can use 
the SN formula to generate predictions for any n-point function of ( to any required 
accuracy compatible with the neglect of spatial gradients. 

The manipulations described above are now standard in the literature. However, 
the foregoing discussion makes it completely clear that the 6N formula is merely an 
identity which expresses ( in terms of geometrical quantities. From this point of view, 
( is nothing more than a book-keeping object which describes how much one region 
of the universe has expanded (if C > 0) or contracted (if ( < 0) relative to the mean 
expansion. Therefore, as has been observed before, it is a pure gauge mode in exact de 
Sitter space [6T] . 

How are we to assign an interpretation to the notion of mean expansion? This is 
necessary to determine what question we actually answer when we compute properties 
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of ( using the 6N formula. It corresponds to specifying which region of the universe 
is supposed to be described by the unperturbed background metric ([1]). In practice, 
we do not imagine that the inflationary patch which gave rise to the large flat region 
we observe today fills our entire spatial shce. For example, this may happen because 
the universe emerged from some primeval phase with initial conditions which were 
chaotically varying, and the requirements for successful inflation were only realised at a 
fraction of spatial positions. Arguments of this sort are invoked for chaotic inflationary 
models, and arise in a similar way for cosmologies based on the so-called "landscape" of 
string theory vacua. In such scenarios we should expect spacetime to retain a disordered, 
inhomogeneous and anisotropic character beyond the boundary of the inflationary patch. 

It is quite clear that we should only associate the FRW background metric ([1]) with 
the region of spacetime containing the primordial inflationary patch, perhaps after a 
few e-foldings of inflation so that the patch is sufficiently close to spatial flatness and 
isotropy [68]. One then expects that whatever primordial unpleasantness lurks beyond 
the boundary is irrelevant when we make predictions from inflation, provided we live 
sufficiently far from the boundary. This is analogous to the decoupling theorem in field 
theory (see, eg., [69]). It follows from the simple point of principle, related to cluster 
decomposition, that in order to make a prediction at one point we should not need to 
know the disposition of the universe at some far distant location. 

In practice this means that a{t) should be interpreted as providing only an effective 
description within some box which coincides with the inflationary patch. But we are 
quite free, if we wish, to imagine that we are working within some smaller box, as long 
as it is large enough to contain scales of order the present horizon size. Whichever box 
size we pick, ( represents the modulation in local expansion as one varies one's location 
within the box. Do the answers to calculations depend on the box size? Certainly. In 
fact, a{t) will generally shift in value between differently sized boxes [TOl [TH |72l [73] . 
because the split between background and perturbations is different on disparate scales. 
One can obtain the correct variation in a{t) by demanding that the expectation value of 
C vanish within the box, that is, ^e{C) = 0, where represents the spatial expectation 
value within a box of characteristic size i. This process, which determines the physical 
perturbation mode within the box, is analogous to imposing a renormalization condition 
which determines a physical property of a particle, such as its mass. 

This ^-dependence of a{t) and ( imply that one is left with predictions for the 
properties of ( which depend on the box size, i. As emphasized by Boubekeur & Lyth, 
one can aggregate predictions in boxes of size i into larger boxes of size i' ^ i in 
such a way that the dependence of C on £ is cancelled by the variation of background 
quantities [TH [73] . This is analogous to renormalization group invariance in fleld theory. 
Nevertheless, the statistical properties of ( itself depend on i. At flrst sight which seems 
like a failure of locality. The resolution is that if the properties of ( calculated within a 
small box are accessible to experiment, the properties of ( calculated some much larger 
box are not, because they require measurements to be made on scales which are larger 
than our present horizon. In order to make predictions for the CMB, we wish ( to 
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represent the perturbation which is actually present in the radiation when the universe 
becomes transparent. This entails making an optimal choice for i in the same way 
that making a prediction for a scattering experiment at some characteristic momentum 
transfer entails making an optimal choice for the running coupling constant. 

In what follows we compute the loop correction in ((^(ki)(^(k2)). This will enable 
us to discuss the ^-dependence of the prediction, and select an appropriate box size for 
the purposes of CMB predictions. 

3. The one-loop correction to the power spectrum of ( 

3.1. A renormalized 6N formula 

In this section, we compute the formal 5N expression for the one-loop correction to 
(C(ki)^(k2)). This expression will turn out to be divergent and will require, as outlined 
in ^ above. In order to organize the calculation, it is convenient to express the various 
contributions to ( using a set of graphical rules, analogous to the Feynman diagrams 
of quantum field theory. Similar diagrams were introduced by Crocce & Scoccimarro 
and adapted for use in the present context by Zaballa, Rodriguez & Lyth [75] . 
The use of such diagrams has recently been formalized to all orders by Byrnes et al. 
[76] using a slightly different notation. Their principal use is to keep track of various 
loop corrections which arise from the 5N algorithm described above; when Eq. (iTOll 
is translated to Fourier space, the products 5(f){yi)^ at a point x become convolutions 
with independent integrals over n — 1 momenta. Some of these momenta are removed 
after taking expectation values, but this cancellation is not always exact and any 
unconstrained integrals can be thought of as making loop corrections to the expectation 
value. 

The rules for drawing diagrams can be summarised as follows. In order to evaluate 
the contribution to an n-point correlator, one draws all diagrams with n external legs 
which inject momenta {kj} into the diagram, where z = 1, . . . ,n. To each external leg 
one attaches a single m-valent vertex (where m ^ 2), at which momentum conservation 
applies. Each m-valent vertex contributes a "coupling constant" d"^~'^N/d(j)'^~^ . The 
internal lines attached to such vertices terminate at any number of special vertices, 
which we label with a cross 'x'. One then draws all combinations of j-gons which can 
be used to connect the x -vertices with no x -vertex left over, replacing each j-gon by 
a j-point expectation value of the 50 as shown in Fig. [H These expectation values 
should be evaluated at the momenta which flow into the corresponding x vertices. For 
clarity, we allow a special notation for 2-gons. (See Fig. [2l) Finally one integrates over 
all internal momenta in order to obtain the expectation value. In such diagrams, it is 
sometimes useful to imagine the various j-gons as representing initial conditions which 
flow from the interior of the graph to the outside, encountering time evolution at each 
m-valent vertex. 

The discussion in ^ implies that we should generalize the formal 6N expansion. 
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Figure 1. Diagrammatic representation of the SN expansion. To represent an n-point 
correlator one draws all m-valent vertices attached to n external legs. Internal lines 
terminate at special vertices marked by a cross 'x '. One then draws all possible j-gons 
which connect the x -vertices, of which two possible examples are displayed. 2-gons 
are treated specially for notational convenience (see Fig. Each m-valent vertex 
contributes a "coupling constant" d™'~^N/d(j)^~^. 

X X h-). 1 I 



Figure 2. Special representation for 2-gon diagrams. In the left-hand diagram the two 
X -vertices should be connected by a 2-gon, but this is difficult to draw. Instead, we 
adopt the notation shown in the right-hand diagram, where the 2-gon is represented by 
two parallel lines transverse to the propagators in which they are inserted. In principle 
there should also be 1-gon initial conditions, where a single x -vertex is tied off at a 
one-point function (Scj)) , but because we are neglecting the one-point function on large 
scales such 1-gons do not appear in these diagrams. 



Eq. ( flOl) . by inserting a zero- momentum expectation value SNq 

C = <5iVo + ^{50.(x)} i^—j N{t,,t,),, (12) 



which is to be adjusted so that lE^(C) = 0, meaning that ( accurately represents the 
perturbation in expansion over a box of characteristic size i. Enforcing this condition 
shows that 6Nq must obey 



^ 1 d^N 

^'^"^-E^a^^W). (13) 

where we have assumed that ((50(x)) = and that statistical homogeneity forces 
(S^)^) = (50-'(x)) to be independent of x for j ^ 2. This guarantees that 6Nq carries zero 
momentum. Angle brackets (■ ■ ■) denote an expectation value in the quantum vacuum 
at past infinity, as usual, and denotes a spatial average over a comoving spherical 
box of characteristic size i. 

The scale factor is not the only background quantity which will receive a 
renormalization. In general we can expect all parameters which specify the homogeneous 
background model to each receive renormalizationsl^l However, the presence of 6No 

f This is similar to the proposal of Kolb et al. [77] , who suggested that a related renormalization of the 
expansion rate at the present day could be responsible for the presently observed accelerated expansion. 
This proposal was criticized on various grounds in Refs. [78l[79l[80l[8Tl[82l[83 l [84 l [85 l [86l[87l[8^ 
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does not affect n-point expectation values of C for > 1; because it leads to 
disconnected contributions. These are not the contributions we wish to measure from 
CMB observations and therefore can be discarded in the present context. 

3.2. The power spectrum of the primordial curvature perturbation 

The simplest possible graph corresponds to Eq. ( fTTl) . which arises from the product of 
the two linear factors in (JTOl) . 

-^^^ (H^) (^'^(ki)'^0(k2))*, (14) 

where the subscript denotes evaluation on the initial slice. Although Eq. f[T^ 
is clearly a tree diagram in terms of the 6N diagrammatic rules described in §3.H 
it implicitly contains loop corrections which are present in the initial correlator 
((50(ki)5(/)(k2)). This illustrates a general feature which arises in calculating ( correlators 
using 6N] there are two distinct types of loop. The first sort are the quantum-mechanical 
loops (called q- loops in Ref. [I6]) present in the expectation values of the 6(f). The 
second sort are loops which arise from expanding ( correlators using the SN algorithm 
(called c- loops in Ref. [16]). To arrive at a consistent result at any given order in the 
loop expansion, it is necessary to retain loop contributions from both these sources. 
Therefore, we require the q-loop correction in (50(ki)(50(k2))=„. This q-loop has been 
computed by Sloth [5T1 [52] in the limit where is long after the mode with wavenumber 
k ^ ki = k2 exited the horizon, and has recently been given in Ref. [58] in the opposite 
case where is the horizon-crossing time associated with k. We are free to base the SN 
expansion on either of these slices, and a discussion of their merits is deferred to §3.31 
below. 

At whatever time it is evaluated, the q-loop associated with Eq. (fT^ is accompanied 
by a number of c-loops arising from 6N integrations. The first of these can be written 

O- - + J (^^'.(kl.k.-q.q). (15) 

where B is the bispectrum of the field fluctuations, defined in Eq. (j4]), and the subscript 
again denotes evaluation on the initial slice. This term gives a contribution to the loop 
correction which depends on the presence of non-gaussianity in the field perturbations 
6(j). We shall see in §4.21 below that this term can in fact give the largest contribution 
to the loop correction if /nl grows on very large scales, although this does not happen 
in the usual scenario. The presence of a single unconstrained momentum integral, here 
labelled q, indicates that this is a one-loop contribution. The finite time of evaluation, 

elsewhere. The present discussion differs from that given by Kolb et al. because there is no suggestion 
here that renormalization of the acceleration rate would have experimental consequences for observers 
within any given horizon volume; indeed, precisely the opposite is the case. The renormalization is only 
relevant for so-called "meta-observers" who can simultaneously compare rates between many different 
Hubble patches; it is relevant for what we can compute, but not necessarily for what we can observe. 
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at conformal time rj^, provides a natural upper cut-off in the momentum integral. On 
scales characterized by wavenumbers larger than q ~ no perturbations have yet 

been generated; on these scales, the spectrum and all higher correlation functions such 
as B and the trispectrum T are zeroJil The upper limit on the integral therefore depends 
on when we take the initial slice to be constructed, whereas the comoving box size £ 
provides a natural (and time-independent) lower cutoff at g ~ 

The next contribution arises from the product of two tree-level spectra. There are 
two ways such a product can arise. The first corresponds to the diagram 

-O- ^ (2-)^^(ki + k,)^ (^) / ^n(|k, - q|)P.(r), (16) 

whereas the second corresponds to a different configuration 

dN d^N r d^Q 
- (2.)3.(k. + k.)^^ J (IT) 

and the remarks below Eq. (ITB]1 concerning the presence of an upper- and lower-cutoff 
in the momentum integral apply equally to both these loops. 

These diagrams exhaust the possible one-loop contributions if we assume that 
the one-point function of the field fluctuations is zero, (50(k)) = 0. As discussed in 
^ the one-point function O is given to one-loop order by a renormalization-scheme 
dependent number which can be absorbed into the background field 4>{t). This number 
is accompanied by other unknown renormalization-scheme dependent numbers which are 
present in the two- and higher ra-point correlation functions. We will assume that to a 
good approximation, all such unknown quantities can be ignored; since we will not keep 
track of scheme-dependent numbers in the higher n-point functions there is no point 
keeping track of O or a renormalization of the background field 0(t) either. Therefore, 
all contributions from one-point insertions will be discounted in the remainder of this 
paper. 

These expressions agree with those recently given by Byrnes et al. in Ref. [76] . 
3.3. Where should the initial slice be set? 

In the discussion of §3.2l the initial time at which the SN expansion is to be constructed, 
r/^,, was left arbitrary. The initial time plays two roles in the analysis: it determines 
whether the field correlators (50(ki) ■ ■ ■ 6(j)(kn))* are to be evaluated at the approximate 
time when the mode with wavenumber A; ~ A;i = crosses the horizon, or if they are 
to be evaluated at a much later time; and it provides an upper cut-off in integrals over 
c-loop momenta. These two effects must combine in such a way that the curvature 
perturbation is independent of rj^,. 

I If we were computing these objects using the underlying quantum field theory of ( rather than by 
following the SN method, then these wavenumbers would lie in a region of Ty-integration where the 
integrand is strongly suppressed by decaying exponentials, leaving almost no net contribution. 



One-loop corrections to the curvature perturbation from inflation 



13 



The loop correction evaluated long after horizon crossing has been computed by 
Sloth [5T| |52] in terms of the variance (50^) = / dlnqV{q), which is in principle a 
infrared divergent quantity and should be cut off at a momentum scale of order the box 
size, g ~ if "P does not vanish as g — )■ 0. The loop correction is found to diverge 
essentially logarithmically with the time of observation, 



where \kri^\ ^ 1 by assumption, a is a numerical factor of order 0(e), is the tree- level 
power spectrum, which is assumed to be almost constant between horizon crossing and 
the time of observation, and Ci(x) is the real part of the imaginary exponential integral. 



Ci(x) = Re <^ Eii(ix) L where Ei„(z) = / — e"*^ 7 -ivr -In^ for z < 1, (19) 



where 7 is Euler's constant. Eq. ( !T9l) apparently contains two distinct divergences. 

One divergence is implicit in (50^). In exact de Sitter space V{q) is a constant, so 
(50^) appears to contain a logarithmic divergence at zero momentum. This is caused 
a pile-up of modes outside the horizon, arising from gravitational redshifting, and is 
ultimately traceable to the quasi-de Sitter expansion. As modes are redshifted outside 
the horizon they carry perturbations, and the ever-increasing number of modes in this 
infrared phase space leads to large effects. However, the logarithmic divergence itself is 
fictional because the original inflationary patch was presumably not infinite in extent. 
Working in a finite box cuts off the divergence on a momentum scale q ~ 

The second divergence is present in the 77* — )■ limit, and becomes manifest when 
one observes the correlator long after the k-mode in question has passed outside the 
horizon. This divergence comes from the logarithmic term in Ci{—2kri^) . As with the 
pile-up of modes, this divergence is fictional in the sense that no observer within the 
universe can really set 77* = 0. Instead, observations must be made at a finite time and 
this finite time gives a cut-off on 77*. 

Although we do not encounter a genuine infinity in either case, these "divergences" 
can generate large coefficients in the perturbation series. Consider first the divergence 
as ?7=K — > 0, which apparently signals a breakdown of perturbation theory at some time of 
order ~ 1/e e- folds after horizon crossing. (Related breakdowns have been reported 
elsewhere in the literature; see, for example, Refs. [46l [26l ETj.) The presence of this 
term endows the correlator with a time dependence that is not visible at tree-level. 
Therefore one should worry whether this time dependence could spoil conservation of 
C on superhorizon scales. If so this would be cause for serious concern, because (- 
conservation is known to hold to all orders of the 6N formula [56] and therefore cannot 
be spoiled by c-loops, and since no field correlator is ever evaluated at the time ( is 
observed it cannot be spoiled by q-loops among the {S(j)} either. 

In fact, the time dependence in Eq. f|T9|) is merely what is required to maintain 
conservation of (. Consider any two widely separated initial times. When evaluated 
at the same time of observation, the coefficient ^N/^(j)^, evolves between these two 
initial slices because the classical field (j) is rolling down its potential. The correlator 
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((50(ki)(50(k2)) must likewise evolve between the two initial slices to compensate, in order 
that ( remains constant. It is this evolution which is represented by the logarithmic 
"divergence" in ( fT9l) . Indeed, Eq. ( lT9i) should be thought of as the beginning of a Taylor 
series in ln(— fc?7*), with the higher-order terms in the series generated by higher-order 
loop corrections. This series is of the form studied by Gong & Stewart [Ml |90]. In 
particular, it implies that perturbation theory in powers of slow-roll parameters breaks 
down at sufficiently late times, of order N ^ 1/e e- folds after horizon crossing. 

Despite this breakdown in perturbation theory, our ability to make predictions 
from theories of the early universe is not really impaired. In such situations, one would 
ordinarily attempt to resum the large logarithms using a form of the renormalization 
group equation. This resummation procedure effectively moves the large coefficients 
from higher-order terms in the perturbation series into the leading-order term. In 
the context of Eq. ( fT9l) . resummation of the large logarithms clearly yields the 
superhorizon evolution of the correlator. Since the perturbations behave classically 
in the superhorizon limit, we should expect that this evolution will correspond to the 
classical evolution of the field [Ml EH |92l |93] . 

Now consider the alternative choice, where the loop correction associated with a 
mode of wavenumber k is evaluated at the time when k crosses the horizon. At horizon 
crossing it follows that —kfj^, ^ 1. The large logarithmic term Ci(— 2A;?7*), and any 
higher-order logarithms, are therefore completely negligible. This version of the loop 
correction has recently been computed in Ref. ^8], with the result 

(<50(ki)50(k2)), = (27r)35(ki + k2)PM (^l-^V.\nk + --^. (20) 

The scale which accompanies k inside the logarithm is set by the details of physics in 
the ultraviolet |§| Although such corrections are extremely interesting, they are not the 
focus of the present paper. For this reason we will ignore them in what follows. At this 
order, (50(ki)50(k2))* is free of infrared divergences. 

The key lesson to draw from a comparison of Eqs. ( fT9l) and ( l20l) is that the 
divergence at late times, associated with powers of ln{—kri^), has been resummed into the 
SN coefficient assocated wth this correlafou fuucfouglt is a remarkable virtue of the 
SN formulation that this resummation is implicit in our freedom to fix the initial slice at 
any convenient time, without any necessity to deploy a complex formalism based on the 
renormalization group. Eqs. ( fT9l) and ( l20l) are therefore merely different representations 

§ In vl-v2 of the arXiv version of this paper, and the version which subsequently appeared in JCAP, 
the regulating scale which accompanies the logarithm was erroneously identified as the infrared scale 
£. I would like to thank E. Dimastrogiovanni for drawing my attention to this error. In these versions 
the coefficient which accompanies In k was also given incorrectly, owing to a calculational error which 
was present in vl-v2 of Ref. [581. This error was corrected from v3 onwards of that reference, to which 
the reader should refer for details. 

II There is nothing to prevent powers of ln(— fc?7*) accompanying q-loop corrections such as Eq. (|20[) . 
These are genuine quantum corrections to the time evolution, and we should not expect the classical 
SN formula to account for their resummation. However, since these terms are suppressed by powers of 
we assume that their contribution is small. 
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of the same physics. It does not matter which time of evaluation we pick, as long as we 
have the ability to perform the relevant calculations consistently. 

In this paper, we take the view that for analytic calculations it is better to evaluate 
the correlators of the fields soon after horizon crossing, where the logarithmic term 
Ci{—2kri^) , and any higher-order logarithms, are negligible. There are two key reasons 
to adopt this approach: 

(i) The q-loop correlators given in Eqs. (fT9|) and (120|) were derived on the assumption 
that slow-roll is a good approximation between the time of horizon crossing and 
the time of observation of the 6(j). At any given order in the loop expansion one 
encounters a potentially large number of diagrams, but many of these diagrams will 
be subdominant when slow-roll applies. For example, as discussed in Ref. [58], the 
q-loop correction to the two-point function has a diagrammatic expansion of the 
form 



loop correction D 




Both of these diagrams are one-loop and therefore of second order in cosmological 
perturbation theory, which is an expansion in powers of 10^^°. However, the 
second diagram is suppressed by an extra factor of order ~ e, since each three- 
point vertex carries a factor of 4>/H [91] whereas the four-point vertex carries 
no suppression [95]. On the other hand, large logarithms such as Ci{—2kri^,) can 
compensate for the smallness of e. Where large logarithms are present the first 
diagram is dominated by the terms of subleading order in slow-roll, since it is 
these terms which are accompanied by large logarithmic factors [511 [52]. When 
terms such as Ci{—2kr]^,) are large one must be wary that the second diagram 
(among other possible sources of subleading slow-roll terms) does not also begin 
to contribute. As the logarithms become increasingly large, further diagrams at 
two- loops and beyond may also become relevant. Therefore, in order to compute 
reliably with a slow-roll truncated q-loop such as Eq. ( 120|) one should pick the point 
of evaluation in such a way that large logarithmic terms do not appear. This means 
that the initial hypersurface on which we base the 6N expansion should be close 
to the horizon crossing hypersurface. The slow-roll expansion is then safe between 
these two hypersurfaces, and we can compute the q-loop to very good accuracy by 
assuming it is dominated by the leading slow-roll part of the first diagram. For 
the two-point function this difficulty is not severe, since for moderate values of 
ln{—kri^) one can account for the contribution of the second diagram, together with 
the subleading part of the first. However, for higher-order correlation functions this 
may become a practical consideration. 

(ii) Eqs. ( !T9|) and fl20l) are computed from integrals over t] whose integrands are rapidly 
oscillating when \kr]\ ^ 1. This kills any contribution from very early times, leaving 
the integral dominated by times where {k-ql < 1. The pre-factors in Eqs. (fT9l) 
and (l20l) are derived by approximating slowly varying quantities such as H and e 
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at the time of horizon crossing. This will be a good approximation if the correlator 
is observed just a few e-folds after horizon crossing but may be on the verge of 
breaking down if In \ kr]^\ is large, since in general one would then expect H and r] 
to evolve between horizon crossing and the time of observation. One can account 
for this evolution using the 6N expansion. 

In the remainder of this paper only the "early" initial slice is considered, where the 
correlator is evaluated just a few e-folds after horizon crossing. 

4. One-loop renormalization of the ( power spectrum 
4.I. The loop correction 

We are now in a position to apply these general principles to the object of central 
interest, the power spectrum (^(ki)^(k2)). After collecting terms in Eqs. (|T4|) - ([T7|) . the 
connected part of the power spectrum can be written 

(C(kOC(k2)) = (27rmki + k2)|(|^yA(A:) + / (§3S(k,q)| (21) 
where E is defined by 

In this formula k can be taken to be either ki or k2, since its orientation is immaterial, 
and the integration d^g is restricted to momenta between g ~ and q ~ £~^. As 
discussed above, the time rj^, at which we construct the initial slice for the 6N expansion 
is taken to be the horizon-crossing time for k. 

To estimate the magnitude of the loop correction in (121]) we must be able to evaluate 
the q integral, but this is not trivial. In general, B{ki, k2, k^) is a complicated function 
of its arguments [9ll [96] and the integral is difficult to evaluate exactly. We can obtain 
a reasonable estimate by applying the method of Boubekeur & Lyth [71]. In this 
approximation, one supposes that the integral is dominated by its infrared singularities. 
Summing over all singularities, one finds that 

PcW = P.w{(0 + n|. (23) 

where is the tree-level power spectrum evaluated at horizon crossing and 11 is defined 
by 

in which, as mentioned above, we have dropped the ultraviolet loop correction of 
Eq. (l20l) . The functions A and B are defined by 

A ^ {5<P') = J -0^,PM and B ^ {h^S<P') = J -^^J^Mq)- (25) 
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The central issue in the remainder of this section, which will also be important for the 
background renormalization described in §4.2l below. is how one should obtain estimates 
of A and B. 

4.2. The approximation of constant tilt 

First assume that the power spectrum and /nl can be approximated by Eq. (iTl) with 
constant tilt parameters n — 1 and nNL- The A and B functions now satisfy 

A{ke)=V,yn-iiki) and B{M) = V.yn^^+n-iik£). (26) 

where i/rix) is Lyth's y-function [72], defined by 

yr{x) ^ I {I - \xr) . (27) 

The parameters /nl, ""-nl and n — 1 are evaluated at k. Eq. fl26l) is approximately 
valid provided that the box size is not too much larger than the lengthscale associated 
with the mode of interest, k~^, so that the approximation of fixed tilt is reasonable 
throughout the box. 

The loop correction is of second order in cosmological perturbation theory, which is 
an expansion in powers of V. However, the coefficients in Eq. ( l24l) generated by the SN 
expansion enter at different orders in the slow-roll approximation. In particular, (A^')^ 
is of order whereas {N"y + N'N'" is of order 1 and f^^hN'N" is of order e^^^. Which 
of these terms is dominant depends critically on the tilts n — 1 and tinl- 

(i) Positive tilt. This is the case where the spectrum decreases on large scales, also 
known as a blue spectrum. It is not generically produced by single-field models of 
inflation. For positive r, the y-function defined in Eq. f l27p satisfies yr{x) ~ 1/r. 
If the spectral and bispectral tilts n — 1 and unl are both positive, then the loop 
correction is very small. 

(ii) No tilt. For r = 0, the y-function behaves like yr{x) ~ In If n — 1 = ??,nl = 0, it 
follows that all terms in 11 behave like logarithms of ki. To estimate its magnitude, 
suppose that the loop correction constitutes no more than some fraction f% of the 
tree-level. It follows that 

n = V,{AN - N,) ^ V^AN < (28) 

where AA^ is the total number of e-foldings between the scale i and the end of 
inflation, and A^^, ~ 60 measures when the presently observable universe left the 
horizon in terms of e-folds before inflation ended. The second approximate equality 
is valid if the loop correction is large, since this entails AA^ ^ A^*. Eq. (128|) 
can be interpreted as a bound on AA^ in order that the loop correction does not 
signiflcantly modify the tree-level prediction. This gives AA^ < (10^/~^P^,)~^, or 
AA^ < 10« for / ~ 1%. 

(iii) Negative tilt. This is the case where the spectrum increases on large scales, also 
known as a red spectrum. It is the generic type of spectrum produced by single-fleld 
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models of inflation. For r < 0, the y-function has asymptotic form yr{x) ~ xl^l/|r| 
when X ^ 1. In this regime, which term dominates depends whether ^ 0: 
the term involving /nl will be dominant if /inl is negative, meaning that the 
non-gaussianity grows on large scalesj^l For the purposes of making an estimate, 
however, it does not matter which term is the dominant one since in either case 
n ~ {kiy, where t > is the dominant tilt, and the difference in pre- factor is not 
significant. It follows that 

n ~ jP^e""', (29) 

where a is a numerical constant which absorbs the pre-factor. Eq. f l29|) gives an 
extremely powerful constraint on how strongly V and /nl can tilt, but unfortunately 
it is untrustworthy precisely in the interesting limit where ki ^ 1 or if the tilt is 
significant. As discussed in ^ in this limit the approximation of constant tilt 
overpredicts the spectrum. In reality the spectrum does not grow like a constant 
power, but only a power of a logarithm. 

Therefore we conclude that the significance of the loop correction depends on the number 
of e-folds of inflation between the fundamental scale, i, and the scale of interest, k~^. 
Except where there is a significant negative tilt the loop correction is small, perhaps of 
order 1% in a model with 10^ e-folds of inflation. However, to study the case of negative 
tilt a more sophisticated analysis is required. 



4-3. Monomial potentials. 

To obtain a better prediction for ki ^ 1, or where there is negative tilt, consider any 
local field theory model of inflation, with potential given by the monomial |97] 

V{<P) = gct>' (30) 

where g is a. coupling constant and r is usually taken to be an integer. For r = 2 and 
small values of the field this corresponds to an Nflation-type scenario [9BJ, in which 
the quadratic behaviour arises as an approximation to a periodic axion potential. For 
larger r or large fields such monomials fit into the general framework of so-called chaotic 
inflation, where is taken to be descending from above the Planck scale with a spectrum 
of field values which may cause some regions of the universe to inflate at the expense of 
other regions in which the field value is not suitable for inflation to occur. 

Inflation occurs for field values greater than 4>f = t/\/2- However, one cannot 
contemplate arbitrarily large values for 0. Local field theory presumably breaks down 
as an approximate description of nature at least when V ^ 1 {m Planck units), which 
occurs at a field value 0p = g^^^^. This is certainly the largest value of the field for which 
the tools used in the present analysis make sense, but there may be other effects which 
imply that our approximations break down even before one reaches the Planck scale. For 

f This is not the usual situation, since in single-field slow-roll inflation /nl oc e [61j , and therefore the 
non-gaussianity decays on large scales. At the times when the fluctuation on large scales was imprinted, 
e was much closer to zero. 
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example, the 'self- reproduction' phase of eternal inflation ends when 250^ > 2e/3 
which implies a maximum field value (p^v ~ {27i'^r'^ g~^y^^^~^'^\ One can imagine either of 
these field values as natural energy scales at which to specify the masses and coupling 
constants of an effective field theory which approximates some ultraviolet completion of 
the Standard Model plus gravity. 

Let us restrict attention to a regime where local field theory is a good 
approximation. Any two arbitrary field values (pa and 06 separated by AA^ e-folds 
of infiation are related via the rule 

cPl = <Pl + 2rAN. (31) 

In order to match the fiuctuation amplitude V( = (24.1 ±1.3) x 10^^i|| observed in CMB 
experiments, the coupling must be tuned to precisely satisfy 

9^^. (32) 

where 0* is the field value when presently observable scales were leaving the horizon. 
This typically occurs roughly N^, ^ 60 e-folds before the end of infiation in order that 
the fiatness and horizon problems of classical cosmology are resolved. Eqs. ( !3T|) -( !32l) 
therefore imply a very small coupling constant, which is hard to understand at the 
Planck scale on the basis of our present knowledge of physics. 

Suppose that infiation begins at some field value 0uv- After a few e- foldings of 
infiation, we can suppose that the fiat FRW metric ([1]) applies within the infiating 
volume. One therefore sets the box size i to correspond to the horizon size at this time. 
The field value when any smaller scale q exits the horizon is related to the field value at 
the onset of infiation by ( I3T1) . giving 

= 0uv-2rlng£. (33) 

It follows that (50^) can be written [99l [51] 

where the second approximate equality applies in the limit 0^, <^ 0uv- If infiation lasts 
for a large total number of e-folds N ^ r, then Eqs. ( pTl) and ( !34l) imply 

This scaling of (50^) = A with the total number of e-folds, namely A oc N'^+^^'^, was 
given earlier by Sloth [52]. A similar procedure allows us to estimate B, giving 

I The quoted value is for WMAP3+SDSS, ignoring any possible tensor contribution POI . 
§ In fact, 0uv should be reduced by a small amount since we are assuming that £ does not correspond 
precisely to the inflationary patch when (p = ^uv but to some slightly later time. However, if the field 
is in slow-roll at this time then (p will not roll far during the few e-folds which elapse. 
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SIGNIFICANCE 


r 

maximum N 


= 2 model 




r : 

maximum 


= 4 model 

iVp 




0.1% 


4.06 X 10^ 






8.65 X 103 






1% 


8.74 X 10'' 


1.28 X 10^° 


5.03 X 10^ 


1.54 X lO'^ 


6.31 X 10^ 


2.50 X 10"* 


10% 


1.88 X 10^ 






2.74 X lO'^ 







Table 1. Limits on the total number of e-folds, N, in monomial models with r = 2 
and r — A. Np measures the number of e-folds between the Planck scale and the time 
presently observable scales were leaving the horizon, and Ngr is a similar measure, 
with the Planck scale replaced by the 'self-reproduction' scale of eternal inflation, (psr 
(defined in the text). 



In making this estimate, we have approximated /nl ~ e at each scale. Although this is 
not strictly numerically accurate [61], it should be sufficient to capture the tilt of the 
non-gaussian fraction. Indeed, Eqs. (I35l)-(l36l) show that A grows faster than _B by a 
power of A^, and therefore that the non-gaussian component is negligible for a large total 
number of e-folds. This is in agreement with our expectation that the contribution to the 
loop correction from any intrinsic non-gaussianity among the 6(j) should be suppressed 
if /nl decays on large scales. In addition, since both A and B only grow as powers of 
N, these equations justify the statements made above that a constant tilt overpredicts 
(50^) for a large range of e-folds. In what follows we assume that N is sufficiently large 
that it is a good approximation to ignore the intrinsically non-gaussian contribution. 

Let us estimate how many e-folds can elapse before the loop correction becomes 
significant enough to spoil the tree-level prediction. The coefficient of the A-term in 
Eq. (jaD is 0(1) in slow-rolljjl 

Thus, in order that the loop contributes at no more than some fraction / of the tree-level, 
we must require Il{dN/d(f)^)^'^ < /, or 

This gives constraints on which are summarised in Table [H Note that the bounds 
arising from the loop contribution are roughly compatible with the bounds derived by 
Wu et al. [26] on the basis that amplification of the quantum geometry should not spoil 
scale-invariance of the spectrum. 

The numerical limits on A^ cannot be taken literally for loop:tree ratios as large 
as 10%, because the estimates for g and (50^) [given in Eqs. ( I32l) and (El])] were 
derived assuming that the tree-level prediction dominated the power spectrum, and 
this is no longer true when the loop constitutes an appreciable fraction of the tree- 
level. Likewise, as the significance of the loop increases, contributions from two-loop 

II Note that in a monomial model there is a special simplification, since N'" — 0. 
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corrections and beyond may become relevant. Nevertheless, these limits show clearly 
that if inflation is taken to persist for an exponentially long time, corresponding to 
initial conditions at the Planck scale or the 'self-reproduction' scale of eternal inflation, 
then the loop correction will be large and cannot be ignored. On the other hand, if 
inflation is only taken to consist of sufficient e-folds to solve the horizon and flatness 
problems, ~ A^* ^ 60, then the loop correction will be completely negligible. These 
conclusions are in qualitative agreement with the results of Sloth |52], who found 
corrections of 0(1) for inflation beginning at the self- reproduction scale. The qualitative 
disagreement between Refs. [511 [52] and the present paper on the one hand and the 
conclusions of Martineau & Brandenberger [87] on the other arises because in Ref. [87] 
the random distribution of phases present in the inflationary perturbations caused an 
overall cancellation of the effect. In the present formalism the distribution of phases is 
invisible; only the spectrum is relevant, not the modes themselves, and the large effect 
described in Table [1] is effectively the phase-space volume divergence given in Eq. (46) 
of Martineau & Brandenberger [87] . 

5. Discussion 

In this paper, the SN formalism has been used to estimate the magnitude of the one-loop 
correction to the inflationary power spectrum in the case of a single-field, slow-roll model 
of the infiationary epoch. The calculation can be split into two separate parts. The first 
depends on a detailed estimate of the one-loop correction generated by interference 
among quantized field modes around the time of horizon crossing. This calculation 
has recently been performed in Ref. [58]. The second measures the pile- up of field 
modes outside the horizon which are generated by the quasi-de Sitter expansion. This 
accumulation of modes as infiation proceeds gives an ever-larger phase space for infrared 
divergences [87], which is measured by the so-called c- loops of the 6N formalism. In 
this paper these two parts are assembled, allowing an estimate of the magnitude of the 
loop correction which includes all relevant effects. 

The results is entirely consistent with our expectations from flat space quantum 
fleld theory: we are free to formulate the fundamental microscopic theory which governs 
inflation at any energy scale we wish, but if we choose to formulate it at an energy scale 
which is wildly different from that of the phenomena we wish to observe, then we must 
expect signiflcant corrections to be generated by loop effects. Such corrections take the 
form of large logarithms that violate the predictions of naive dimensional analysis. In 
the present case the loop correction is given by Eqs. fl23l) - fl2^ . the fundmental theory is 
formulated at a scale i, and we wish to predict phenomena associated with the CMB at 
a scale k~^. The large logarithm which gives a signiflcant loop correction is iV ~ \nk£. 
If ki is exponentially large, which occurs in monomial models of inflation if initial 
conditions are set at around the Planck scale, then the loop correction has a sensitivity 
to the tilt of the fleld perturbation power spectrum. The 6N loop from the convolution 
of two power spectra is the most sensitive and scales like A^^+^Z^; the 5N loop from any 
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intrinsic non-gaussianity among the 6(j) is suppressed, because /nl scales in the opposite 
direction, giving case of a negative tilt, which produces perturbations of 

increasing amplitude on large scales, the loop correction rapidly becomes very significant. 
The correction is irrelevant only if inflation lasts for a small number of e-folds. 

The physical effect which demands that loop corrections be taken into account is 
nothing more than back reaction. Long ago, Salopek &; Bond |100] used the stochastic 
formalism to show that any initially homogeneous background field would develop wild 
fluctuations on ultra-large scales as inflation progressed. The large loop correction 
encountered in the limit of a large box is another manifestation of this effect. It 
measures the dispersion generated by comparing a correlator calculated in one region 
of the universe with the same correlator calculated in a different region, perhaps where 
the background fields have very different values. The correlator calculated in the large 
box is therefore not measureable in any experiment: it involves comparing the CMB 
temperature which would be seen in an ensemble of different universes, and no physical 
observer can make such a comparison. For this reason, the correlation is a "computable" 
in the sense of Witten [23], but is not an observable. 

This result has no consequences for how we make or report cosmological 
observations, but has very significant consequences for how such observations are 
compared with theoretical models. In models where the fundamental theory is not taken 
to be valid within an exponentially large region, we are entitled to entirely disregard the 
question of loop corrections, since in this case they yield a negligible shift in tree-level 
values. This class of models includes the familiar single-field examples where inflation 
ends when the field reaches a certain point on its potential. In such cases we know 
the conditions which lead to the end of inflation and we always have the option to 
work within a small box|il On the other hand, where isocurvature fields are present, 
these fields will generally acquire large fluctuations within the inflating region. If these 
fields play a role in ending inflation, then we may lose the option to compute within 
a small box, whereas the correlator within a large box may be afflicted with large 
loop corrections. This problem is perhaps most severe precisely for those models to 
which we presently attach the greatest interest - those motivated from fundamental 
particle physics. For example, one popular class of such models arise from the so- 
called landscape of string vacua |101] . In such models, one imagines the unwanted six 
dimensions predicted by string theory to be wrapped up on a Calabi-Yau three-fold 
which may be topologically quite non-trivial, as measured by the existence of a large 
number of homological cycles. These cycles can be wrapped by fluxes of p-form flelds, 
giving rise to regions of the three-fold which have different values for the cosmological 
constant A, together with the other constants of nature |102j . and in general many 
scalar flelds will be present which measure the size and shape of the Calabi-Yau three- 

f Indeed, in such a case, it would presumably make no difference even if we chose to work within the 
large box because the only distinction between widely separated locations would be a time delay. This 
can be gauged away by a reparamctrization of time, meaning that the difference is not physical. This 
is the essentially the argument given by Unruh |35j . 



One-loop corrections to the curvature perturbation from inflation 



23 



fold. These different regions may exhibit substantially different low-energy physics. In 
such models, the natural scale at which one computes the masses and couplings in the 
low-energy theory would be the Planck scale, or perhaps the self-reproduction scale of 
eternal inflation. Eqs. (!23|) -( l2^ show that to compare any such theory to observation, 
one must account for large loop corrections which are inevitably encountered in passing 
from the Planck scale to the scale relevant for CMB predictions. 

This has some implications for how large- field models of inflation are compared with 
CMB observations. To understand how the Planck scale theory is related to the theory 
when scales associated with the CMB exited the horizon, one would need to integrate 
out the relevant degrees of freedom between Mp and the CMB scale. By analogy 
with conventional field theory, as one changes scale the effective lagrangian should be 
controlled by a flow equation analogous to Polchinksi's exact renormalization group 
equation |103] . In fiat space theory this flow has well- understood effects, attracting 
the lagrangian to an invariant hypersurface on which non-renormalizable couplings are 
suppressed. (This would refer to the effective field theory for (, not to the couplings 
in the matter theory.) Unfortunately a complete understanding of this process in the 
cosmological case is not yet in hand, and it does not appear to be understood precisely 
how the Planck scale theory will be related to the theory in a small box. One possible 
approach is to use the method of Salopek & Bond |100j . who gave a Fokker-Planck 
equation which determined the probability distribution of the classical fields on large 
scales. Having done so, it is possible to use a small box to determine the value of the 
correlation functions at each value of the allowed background fields, but one is then 
restricted to making probabilistic statements concerning what should be observed in 
the CMB. This issue deserves further study. 

In any case, considerable work remains to be done. The present analysis should be 
generalized to the case of multiple scalar fields in order that it can be applied to the 
best-motivated models from fundamental physics. This involves an extension of the q- 
loop calculation given in Ref. [58] to the case of many scalar fields. The calculation given 
by Weinberg [H] shows that where ( is accompanied by N" isocurvature fields, we can 
expect the ( loop to receive an enhancement by a factor ~ A/". However, in a multiple 
field scenario a more involved process is required to estimate (50^) |104l \105\ \10Q\ 1107] . 
since the growth of V on superhorizon scales is no longer simple to predict. These issues 
should be addressed with some urgency. 
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